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Abstract 

A classical hydrodynamic approach is used to calculate the magnetoplasma 
excitations of two vertically coupled electron dots. The electrons are con- 
fined by different parabolic potentials in which case Kohn's theorem is no 
longer valid. The equilibrium density profiles of the electrons in both dots 
are calculated as function of the interdot distance. We find that for unequal 
confinements of the two dots the electron density in one of the dots becomes 
ring-like. The electron densities are then used to obtain the magnetoplasma 
frequencies. The oscillator strengths are calculated, and we find that other 
then the center of mass modes can be excited due to the electron-electron 
interactions. 

PACS numbers: 73.20.Dx, 36.40.Gk 
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I. INTRODUCTION 



Modern semiconductor technology makes it possible to fabricate low-dimensional nanos- 
tructures with controllable chemical composition and geometric structure. In artificial atoms 
a number of electrons are confined into a single quasi-two-dimensional quantum dot which 
has been widely studied over the past few years (see, for instance, Ref. |l]). Most of the work 
has focused on quantum dots with a soft lateral confinement potential in the heterostructure 
plane (2D-plane) which are the so called parabolic quantum dots. Such quantum dots ex- 
hibit a far-infrared absorption spectrum which refiects only the center of mass motion of the 
electron system which is a consequence of a generalized Kohn's theorem.! When a magnetic 
field is applied perpendicular to the plane of the quantum dot it results into two absorption 
peaks whose positions are independent of the number of electrons in the dot. Note that 
the same absorption spectrum follows from a pure classical analysis. In order to reveal ef- 
fects due to electron-electron interaction the confinement potential should be anharmonic, 
or band non-parabolicity effects should be taken into account. 

The other possibility to trace the influence of electron-electron interaction on the quan- 
tum dot properties is to consider the more sophisticated system of coupled quantum dots. 
An example of such a system are two vertically coupled dots. The properties of such a sys- 
tem for small number of electrons were studied previously via exact diagonalizationi or for 
a model system with an exactly solvable interaction! where the main attention was on elec- 
tron correlations in the ground state, magic numbers and the possibility of the appearance 
of new type of entangled electronic states. In Ref. |^ it was shown that in the case of a double 
dot system, containing three electrons, with different parabolic conflnement potentials, the 
electron correlation effects exhibit themselves in the absorption energy. The classical double 
dot system containing a flnite number of electrons exhibits a rich variety of flrst, second and 
higher order transitionsJl Two vertically aligned dots is the flnite system analog of the well 
studied bilayer systemJ^ which have been shown to exhibit much richer physical phenomena 
then 2D and 3D systems. 

Such coupled dot systems can be realized for example in vertically etched double quantum 
well structures,! parallel layers of self assembled quantum dotsi and by using a conducting 
atomic force microscope tip above a quantum well which can induce locally a dot structure 
in the quantum well and one beneath the lowest barrier.0 Recently, the successful growth 
of vertically aliened dots of up to 10 layers of InAs islands separated by GaAs spacer layers 
were reported.til These experimental attempts are still in their infancy and to our knowledge 
no magneto- absorption experiments have been performed yet. 

The purpose of the present paper is to study the properties of the system consisting 
of two vertically coupled quantum dots in a uniform perpendicular magnetic fleld. We 
calculate the electron density proflles, the spectrum of magnetoplasma excitations and the 
corresponding power absorption. We restrict our consideration to dots with a large number 
of electrons. In this case the classical hydrodynamic approach is adequate, and we follow the 
ideas presented in the paper of Ye and ZarembaEl where the single quantum dot with quartic 
conflning potential corrections was considered. A similar technique was used to study one 
layer electron dots with parabolic0'Ei and non-circulaJii conflnement potentials. 
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II. MODEL 



We consider two infinitely thin two-dimensional circular symmetric dots which are ver- 
tically aligned and are separated by a distance a as is schematically shown in Fig. |l|. The 
electrons are confined by different parabolic potentials. In the hydrodynamic approach the 
electrons are described by density functions Pj(r) where j = 1,2 refers to the corresponding 
dot, and r = [x, y) is the position coordinate in the plane of the dot. Those electron densities 
create an electrostatic potential (p{Y, z) which obeys Poisson's equation 

VV = —{pi{v)5{z) + p,{v)5{z - a)}, (2.1) 

where e is the dielectric constant of the medium the electrons are moving in. We assume 
that e is constant thrroughout the whole system which is a good approximation for e.g. the 
GaAs/Al^.Gai_a;As system. 

In a purely classical approach, the electrons in both dots achieve an equilibrium distri- 
bution in which the lateral force arising from the confining potentials is balanced by the 
electrostatic repulsive force due to the electrons. Equivalently, the equilibrium conditions 
can be expressed as the constancy of the chemical potential in each of the dots given by the 
expressions 

ki 

ii\ = —r"^ - eip{Y,z%=Q, (2.2a) 
P2 = —r'^ - eip{Y,z%=a- (2.2b) 
Here the strength of the parabolic confinement potential is characterized by the correspond- 



ing factor ki which is related to the confinement frequency as uoj = ^Jkj/m. Here and further 
the symbol m stands for the effective electron mass. The difference in the above chemical 
potentials can be controlled by applying some external bias potential 

eV = pi-p2- (2.3) 

In equilibrium the electron densities are circular symmetric Pj(r) = Pj{r) and vanish 
outside some radius Rj which must be calculated self-consistently. Those radii are related 
to the number of electrons in the dots 

2ti / pj{r)rdr = Nj. (2.4) 



Condition ( p. 3D together with definitions (|2.2| ) and Poisson's equation ( p.l|) form the basic 
set of equations from which the equilibrium electron density profiles in the dots can be 
obtained. In a classical description the equilibrium electron distributions are not infiuenced 
by an external magnetic field. 

In order to obtain the excitation spectrum and the power absorption in the coupled dot 
system the small harmonic deviations from the electron equilibrium density and potential 
functions have to be introduced 

5p,(r), 6^ir)^e'-\ (2.5) 
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Those additional quantities obey a set of dynamic equations. They are composed of the 
same Poisson equation for the above deviations and the continuity equations 

iu;e6pj{r) + V ■ij{r) = 0. (2.6) 

Usually the currents in the dots jj(r) are defined by solving the linearized hydrodynamic 
equations for the averaged electron velocities. The solution of these equations leads to Ohm's 
law 

},{r) = -a'" ■V6if{r,z)U=o,a. (2.7) 

In the presence of an external applied perpendicular magnetic field, the components of the 
conductivity tensor cr*^''^ in polar coordinates are given by 

^'■^ - - mico^-uj^J ^^-^^ 

Here the symbol Uc = eB/mc stands for the electron cyclotron frequency, and u is the 
frequency of the external perturbation, i.e. the far-infrared radiation. Poisson's equation 
together with Eqs. ( p^ ) and ( |2.7D compose the basic set for the coupled dots, from which 
the excitation spectrum will be obtained. 

In order to calculate the power absorption spectrum the same equation set has to be 
used but with the single replacement Sep — > 6ip + (p'^^'^, where the external field potential <^^^* 
is added. Then the frequency dependent power absorption spectrum is calculated in the 
standard way using the following expression 

P{u) = -- dt Y.J d'ri, . V^^^*. (2.9) 



III. EQUILIBRIUM PROPERTIES 



The first step in the static problem is obtained by replacing Poisson's equation (|2.1| ) by 
the homogeneous Laplace equation which is solved with the following boundary conditions 
for the potential in the dots 



dip 
dz 



27re 



-Pi, 



z=0 



dip 
dz 



27re 



-P2- 



(3.1) 



As was shown in Refs. |r^]-[T^ for the case of a single dot the most simple way to treat 
this static problem is to use the oblate spheroidal coordinate system centered in the dot, 
which is defined as follows: x = ^ (cr^ + 1) (1 — r^) cos 9, y = ^ (a^ + 1)(1 — r^) sin6', z = 
ar. In this coordinate system the disk corresponds to the region a = 0, < r < 1. 
The advantage of this coordinate system is that the variables in Laplace's equation can 
be separated and the solutions can be expressed as a linear superposition of the potential 
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harmonics0 p\p\{^r)Qf\ia), with J > \m\. /^''"I and Q^^ are the associated Legendre 
polynomials of the first and the second kind, respectively. Now the static problem is reduced 
to the solution of an algebraic set of equations for the expansion coefficients which follows 
from the equilibrium conditions. 

We shall use a similar approach in our coupled dot problem. In order to have the 
most simple expressions for the main part of the electron-electron interaction, namely the 
electron-electron interaction in the same dot, we introduce two coordinate systems each 
centered in the corresponding dot. We rescale both coordinate systems, replacing r — ^ ri?,, 
and introduce the oblate spheroidal coordinates in each of the dots, respectively. Next, we 
divide the potential into two parts ip = (pi + ip2, each created by the electron density in 
the corresponding dot, and introduce dimensionless variables by performing the following 
transformations 

Pj{r) — > pj{r)pfVl-r^, <Pj — > ^\^^^- ^^''^^ 

By this transformation we explicitly took out the non-analyticity of the density function 
at the dot edge which is inherent in the case of a soft confining potential.li3 If we choose 
the scaling factors p^p = 2ekjRj/'rr'^e'^, introduce the dimensionless numbers of electrons in 

the dots rij = 3Nj/27rR'^pf\ chemical potentials vj = pj/kjB?-, and ratios 7 = R2/R1 and 
K = k2/ki, we transform our static Eqs. (|2.2| ), ( |3.1|) , and ( p.4|) into the following set of 
dimensionless equations 



da 



1 

2r2 



r'/2 + <*(r) + YK^T{f) = vu (3.3a) 



rV2 + —vT\r) + y^f (r) = v^, (3.3b) 



(T=U 



Pj(r), (3.3c) 



/ drr^/l — r'^pjir) = / drr'^pjir) = —. (3.3d) 
JO Jo 3 

Here the symbol f stands for the coordinate in one dot calculated in the coordinate system 
of the other dot. The superscripts indicate whether the potential is created by the electron 
density in the same dot ('^°*) or the electron density outside it (°"*). 



Laplace's equation must not be inc 



above, its solutions are already known ilB and can be presented as 



uded into the set of equations because, as is discussed 



^rw = E^F^^2K-)7^, (3.4) 

where due to the circular symmetry of the equilibrium state only the even potential har- 
monics with m = have to be taken into account. This general expression can be used for 
the potential outside the dot. Inside the dot (cr = 0) the potential can be presented in a 
more simple form 

^f\r) = j:cl'^P2i{r). (3.5) 
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The electron density in the dot can also be presented as an expansion in Legendre polyno- 
mials 

P.W = :;^E^?^2z(r). (3.6) 
I 



Here, the additional factor is included in order to simplify the boundary condition ( p.3c| ) 
which now leads to the relation = Lj^C/"''' between the coefficients in both expansions 
3D and ra, with 



TT r((/ + 


\m 


l)/2 + i)r((/-| 


m|)/2 + 1) 


2r((/ + 


m 


+ i)/2))r((/- 


\m 


+ l)/2) 



Lr = o ;a; i", "w;^;;v^, i' , . (3-7) 



where the symbol T stands for the Gamma function. Inserting expansions (|3.4|) , ( p.5|) and 
( p.6[ ) into Eqs.( |3.34 |3.3b| ), multiplying them by -P2z(t), integrating over r and using the 
Legendre polynomial orthogonality condition 

Jo = 2l + l{l-m)V ^^-^^ 

we obtain the following set of linear equations 

Cl'^ + KY.TUlh,dMCi^^ = (v, - 5,,o + ^5/,2, (3.9a) 

Here the symbol T^^ stands for the overlap integral of the potential harmonics created by 
the electron densities in the different dots 

TUl,d) = / dTP2l{T)P2n{r) , 3.10 

7 -^0 Q2n{0) 



withf = { [Vg" + - q] /2y^\ a= { [Vg" + + q] /2}''' and g = 7^(1 -r^) + ^2 _ i. 
In the above expressions we used the dimensionless distance between the two dots d = a/Ri. 
Eqs. ( p.9|) should be supplemented by two more conditions (as was shown in Ref. |12|) 



1/2 



5]L2/P2/(0)CP = 0, (3.11) 

which guarantees the non singular behavior of the density functions at the dot edge r = 0. 
Conditions ( p.ll| ) together with Eqs. ( p.9|) form a set of coupled linear algebraic equations 



for the coefficients vj and C/"'"', which is equivalent to a nonhomogeneous matrix equation for 
the vector X = {t>i, Cq^\ ■ ■ ■ , f2, Cq'\ ■ ■ ■}. Only in the case of well separated dots (a — > oo) 
an analytical solution can be found. In that case we obtain vi = V2 = ni = n2 = 1 
which corresponds to the single parabolic dot solution.^ For the general case we solved 
the problem numerically by means of Q-R-decomposition of the truncated matrix. We 
restricted our consideration to 9 potential harmonics for each dot which was sufficient to 
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obtain an accuracy better that 0.1%. We should note that in the case of different confining 
potentials k ^ 1 or with an applied bias eV ^ the numerical solution can not be found 
straightforwardly because of the unknown ratio 7 of the dot radii. We found 7 by solving the 
boundary condition (|2.3|) iteratively. After the ratio of the dot radii is found all quantities 
of interest can be calculated just using the scaling factors introduced in Sec. ^ namely, 

N, = N, N2 = 7=^«:— iVi, 

Til + y^Kn2 rii 



1 

{rii + 7^Kn2)(-^/^) 



-Rl — 7 ; z wi R2 — iRli 



(0) 1 (0) (0) /o 1r,^ 

Here the symbol N = Ni + N2 stands for the total number of electrons in boths dots, and 
two scaling factors — Rq = (Swe'^N/AekiY^^ and po = 3N/2nRl — are used. The first one 
corresponds to the equilibrium radius of the first dot when all electrons are situated in this 
dot. 

The numerical results illustrating the equilibrium electron densities in the dots are shown 
in Figs. and ^ The density profiles for equivalent coupled dots = 1) are given in Fig. |^ 
for various distances a between the two dots. In the case when there is no external bias 
(yUi = fi2) both dots have the same density of electrons (A^^i = N2 = N/2). In the limiting 
case of small distances between the dots the result reduces to the single parabolic dot with 



density profile pi ~ y 1 ^ v'^/Rq. In the opposite case of large interdot distances we again 
have two parabolic dots where in the above density functions the radius is scaled by 2^/^. 
For intermediate distances, however, the density profile differs from that of a single parabolic 
dot. This is illustrated in Fig. |^ by the thick dashed-dotted curve which shows the parabolic 
dot density with the same dot size and the same total number of electrons as those for the 
coupled dots with interdot distance a/Ro = 0.3. The dependence of the dot size, Ri, on the 
interdot separation is shown in the inset of Fig. With decreasing interdot separation the 
size of the dots increases and the electron density profile is more spread out. In Fig. ^ the 
density profiles for the case of different confining potential strengths (k = 2, pi = P2) are 
given for the same values of the interdot distances as in Fig. |. Notice that now the density 
profiles differ essentially from the single parabolic dot solution. In the case of small interdot 
distances even a ring like distribution can be obtained when the electrons in the dot with 
the weakest confinement are pushed outside. The above ring like distribution effect can be 
even more reinforced by applying an external bias. Notice also that with decreasing interdot 
distance electrons from dot 2 are moved into dot 1 in order to keep the chemical potential 
the same (see inset of Fig. 0(b)). The inset of Fig. ^(a) shows the dependence of the radius 
of the two dots on the distance between the two dots. 



IV. MAGNETOPLASMA EXCITATIONS 

To solve the dynamic equations (|2.1| ), ( p.6|) and (pTTp we use the same oblate spheroidal 
coordinate technique as in the static case which was described in previous section. Restrict- 
ing ourselves only to those modes which can be excited in the dipole approximation we use 
the following expansions for the harmonic deviations of the electron density 
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^P.W = 4E^??AmMe^"'> (4.1) 

I 

the potential inside the dot 

<^^f(r)=EcF^^2Vi(r)e^"^ (4.2) 
I 

and the potential outside the dot 

(r) = E 4'^P2\^iir)^k^e-'. (4.3) 

The symbol a = ±1 is the angular momentum of the excited harmonics actually indicating 
the polarization of the circularly polarized electric field. The above electron density and 
potential harmonics obey the same boundary conditions in the dots as given by Eq. ( p.3c|) 
which leads to the relation g^^^ = Lg^+icf^ We will use the same dimensionless variables as 
in previous static case scaling additionally the frequencies u ^ u ■ ui and ujc ^ ^-^c - '-^i- In 
those dimensionless variables the continuity equation (p]6|) becomes 



Luiu' - ui)T6pi - V ■ ( s^'^ ■ V((5(/?f * + 6<f^') ] = 0, (4.4a) 



^(^2 _ ^2j^^p^ _ V . (^^(2) _ ^ ^ ^^^^^^ 

where for the sake of convenience we introduced the tensor s , which is proportional to the 
static densities, with the following components 

4r = 4e = ^^Pi(^)> s2 = -4? = -iuJcTpiir). (4.5) 

Now inserting expansions ( [4.1| ), ( [4.2| ) and ( ^.3| ) into Eqs. ( |4.44 [4.4b| ), multiplying the 
obtained expressions by rF/(r) exp(— ia^), integrating over 6 and r, and using the orthog- 
onality condition ( |3.8|) we obtain the following set of equations for the coefficients c[*^ 

- ^^^) 2(^ + l)(2^+l) ^i^^^,a) _ ^ _ ci/7)c(f) = 0, (4.6a) 

' ^;) 2('+ ;)P' + l) Li,^,e!" - i i: Hl'(7, d)cW - J: V,M' = O, (4.6b) 

which is cast into a homogeneous matrix equation 

M{uj)X = (}, (4.7) 

for the vector X = {c^\ c[ \ ■ ■ ■ , c[\ ■ ■ ■}. The dynamic matrix A4 is composed of the 
overlap integrals of the potential components belonging to the same dot 

V,^^ = ^ r de C drrP' {r)e-'^'V ■ ( 7^'^ ■ VP' Me'"' 



[drPl,^,{r)\^-'-^[n{n + l)-l] 

PUAr). (4.8) 



-(1-r^ 



dr \ T 



d . ( ds^^l 



dr \ dr 
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and to the different dots 



Wt:^il,d) = t drPl^,{f)^±^ |-^[/(/ + 1) - 1] 




(4.10) 

The functions f, a and q are defined as in the stationary case. 

The magnetoplasma excitation frequencies ujj are obtained as solutions of the following 
equation 

detM{uj)=Q. (4.11) 

The matrix determinant was calculated numerically by means of LU-decompositionlll trun- 
cating the matrix up to 7 potential components for each dot which guaranteed the same 
accuracy as in the equilibrium case. 



V. POWER ABSORPTION 

We restrict our consideration to the case of weak dissipation when the magnetoplasma 
excitation spectrum can be most easily revealed. In that particular case the power absorption 
can be calculated using the same set of dynamic equations ( [4.4a| , [4.4b| ) but now with the 



external electric field added. Our system is circular symmetric and therefore it is most 
convenient to consider the circularly polarized electric field E^-^^ = E^iex + iaey) exp(iujt) 
with the corresponding external potential amplitude (p^'^^{r) = —Eo{x + iay). Now inserting 
the external potential into the set of dynamic equations ( [4.44 |4.4b| ), performing the same 
procedure as before, and scaling the vector X X\/27i\e\EQ/kiRi we arrive at the following 
nonhomogeneous matrix equation 

M{Co)X = B. (5.1) 
Here, the vector B = {&o^'*, • • • , ftg^-*, ■ ■ ■} is composed of the following components 

b? = 2{l + 1){21 + l){u: + au^) Y: S,,ir4> ■ | ^ /j^^ ] Z 3' (5-2) 



and 



Si',i = / —P2i'{r)P2i+i{r) (5.3) 
JO r 

is an overlap integral of even and odd Legendre polynomials. 

Besides, in order to take some small dissipation into account we applied the standard 
procedure making the following frequency replacement u ^ uj = u + it]. The symbol rj 
represents a phenomenological inverse electron relaxation time, and takes into account the 
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internal dissipation within the electron fluid in the dot. According to expression (|2.9|) the 
power absorption can be presented as 

P{u) = -uPolmLlic'i'^ + /€7'c!,'^}, (5.4) 

where we have introduced the absorption power scaling factor Pq = AeEQRluJi/SiT. Formally, 
the dimensionless absorption power can be expressed in matrix form 

p{io) = P{uj)/Po = -Im {T ■ X}, (5.5) 

where the vector T = {ujLI, 0, ■ ■ ■ , 0, UK'y'^L\, 0, ■ ■ ■ , 0} has only two non zero components. 

Equation ( |5.1| ) and the power absorption ( p.5| ) was solved numerically using the SVD 
(singular value decomposition)^! of the dynamical matrix 

M{u) =UA{u)V'^, (5.6) 

where U and V are unitary matrices (superscript stands for transposed matrix) and A is 
the diagonal matrix composed of the eigenvalues Xi{uj) of the dynamical matrix Ai. Due to 
the decomposition ( ^.6|) the solution of equation (|5.1| ) can be presented as 

X = VA~\Co)U'^B. (5.7) 

Inserting it into expression ( p.5|) we obtain the general expression for the power absorption 

p{uj) = -Im {T ■ VA~\(Jo)U^B}. (5.8) 

This expression can be further simplified taking into account that we are considering 
the case of weak dissipation [r] 0) only. In that case the power absorption consists of 
narrow Lorentzians located at the magnetoplasma excitation frequencies uj. According to 
Eq. ( |4.11| ) the corresponding conditions Xi{uJi) = should be fulfilled. This enables us to 
use the following expansion 

Xi{oj) = X[{uJi){uj - Lj, + ir]). (5.9) 
Inserting it into expression (p. 81) gives the following final expression for the power absorption 



p(.) = -Im E ^,^^^1^^ _ ^_ ^ + (5.10) 

where we introduced the coefficients 

/* = Re rjj— , (5.11) 

and the symbol Vi stands for the i-th vector-column of matrix V, and the symbol for 
the i-th vector- row of the transposed matrix U^. Note that we neglected the imaginary 
part of the right hand side of expression (|5.11|) , because for weak dissipation it only leads 
to a negligible eigenfrequency shift. For the same reason the other quantities (like U, V, A') 
should be estimated a.t u = Ui. 

In the case of weak dissipation the quantities ( p.ll|) together with the mode frequencies 
cUj characterise the power absorption. By analogy to quantum mechanical problems they can 
be referred to as the classical oscillator strengths. As is shown in Appendix ^ these oscillator 
strengths obey a general sum rule like their quantum mechanical counterparts, indicating 
that the interaction variation in the system leads to the redistribution of power absorption 
among various vibration modes while the total power absorption remains unchanged. 
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VI. RESULTS AND DISCUSSION 



For the case of equivalent dots k = 1 and /xi = /i2 the eigenfrequencies and oscillator 
strengths are given in Fig. ^. The frequencies in the limiting case oi a ^ oo are shown in 
Fig. ^(a) as function of the applied magnetic field (tUc). The eigenfrequencies for clockwise 
circular polarization {a = 1) are shown by solid curves while the dashed curves indicate 
the eigenfrequencies for the opposite polarization (a = — 1). In the limit case a ^ oo of 
decoupled dots one finds that = and VjJ;^ is a diagonal matrix (see Eqs. (4.9) and 
(4.10)). Then the eigenmodes are the basisfunctions -P2«+i(''") ^^"^ spectrum branches 
can be uniquely labeled with the numbers {I, a = ±1). In this limit the eigenfrequencies 
are twofold degenerate and coincide with the single parabolic dot frequencies presented in 
Ref. |l^. In the low frequency region a large set of branches are found which correspond to 
the edge modes. They are not resolved and shown up as a shaded area. In Fig. ^(b) the 
spectrum for an interdot separation of a/Ro = 0.5 is presented. For the sake of convenience 
we still label the branches by the previous /-values in order to show the correspondance of 
the branch to its counterpart in the a — oo limiting case. In this case in fact several /-values 
are mixed. Note that when the distance between the layers decreases the degeneracy is 
lifted and all branches are split into two and are shifted down in frequency except for the 
two branches depicted by the thick curves which correspond to the center of mass motion 
of the whole system of electrons in both dots. Only the latter two frequencies have non 
zero oscillator strengths which do not depend on the interdot distance and are shown in 
Fig. ^(c). This is in accordance with Kohn's theoremi which can be easily generalized to 
our vertically aligned two dot system for the case of identical dots, i.e. k = 1. The decrease 
of the eigenfrequencies with decreasing interdot distance is a consequence of the increase of 
the repulsive interaction of the electrons in the different dots and is discussed in more detail 
in Appendix |B[ 

The case of nonequivalent dots is illustrated in Figs. ^ and ^ where the results for 
K = 2 and fii = fi2 are presented. The frequencies for the a —>■ oo limit are shown in 
Fig. |^(a) , and the corresponding oscillator strengths are shown in Fig. ^(b). In this case 
the system consists of two noninteracting parabolic dots where the frequencies are scaled 
by the factor a/k = ^72. The branches are labeled with the corresponding /-value. The 
non-zero oscillator strenghts are those which corresponds to the center of mass motions of 
the electron systems in the separate dots. Note that the oscillator strengths obey the sum 
rule that the sum of the oscillator strengths equals 1 independent of Uc- The redistribution 
of the oscillator strengths between the separate dot branches actually reflects the electron 
number redistribution between the dots (see Fig. It is also evident that those spectrum 
branches which are closer to the cyclotron resonance branch {u = Uc) have essentially larger 
oscillator strengths. 

The spectrum branches and the oscillator strengths change fundamentally when 
nonequivalent dots are put closer to each other as is shown in Fig. |^. Now all the branches 
are shifted down in frequency and non- zero oscillator strengths corresponding to other then 
the center of mass motion branches start to appear. For instance, in Fig. ^(b) the oscillator 
strengths corresponding to the frequency branches (/ = 1) are non-zero as shown by the 
curves labeled 5 and 6. The reason is that for nonequivalent dots the electron density distri- 
butions are far away from the parabolic dot distribution (see. Fig. H), and Kohn's theorem 



11 



is no longer valid. Note also from Fig. ^(a) that the frequency shift to lower frequencies is 
much more pronounced as compared to the case of equivalent dots (Fig. §(b)). The above 
frequency shift can be even more amplified when an external bias is applied. This is illus- 
trated in Fig. 1^ where the results for k = 2 and 1x1—^2 = —0.2 are given. From Fig. ^(a) we 
see that now the (/ = 1) branches start to interfere with the (/ = 0) ones. Although there 
is still no intersection with the branches of the same polarization the branch (/ = 1, a = 1) 
approaches the cyclotron resonance line {uj = uJc) so closely that the oscillator strength (see 
Fig. 0(b)) is transferred from the main branch {I = 0, a = 1). 

We would like to note that in the above case of nonequivalent dots when the system no 
longer obeys Kohn's theorem the edge modes have small but non-zero oscillator strengths. 
In order to reveal those weak modes we performed a more accurate dynamic matrix diag- 
onalisation for the case fii = fi2, a/Ro = 0.5 and k = 2 by restricting ourselves to the 
/ = and / = 1 components for each dot. The results for the edge modes (curves 9 and 10) 
together with the oscillator strengths for the upper / = 1 modes (curves 7 and 8) are given 
in Fig. |[ We see that the oscillator strengths become very small for spectrum branches 
which are substantially different from the above mentioned u = Uc. 

VII. CONCLUSIONS 

In this paper, we have studied the magnetoplasma excitations in two vertically coupled 
electron dots in the dipole approximation. We used the classical hydrodynamic approach 
which is adequate in the limit of a large number of electrons. The confinement potential was 
taken parabolic in both dots. We calculated the equilibrium electron density distribution as 
function of the interdot distance. As function of the magnetic field the frequency spectrum 
and the corresponding oscillator strengths are calculated. 

In the case of equivalent dots, i.e. both dots have the same parabolic confinement 
strength, the equilibrium density distribution is similar to the single dot parabolic one, and 
due to the generalized Kohn's theorem only the center of mass motion can be excited. 

If both dots have a different parabolic confinement the electron densities can differ ap- 
preciably from the parabolic case. Even ring like electron distributions are found. Due to 
the fact that both dots have different confinement potentials Kohn's theorem is no longer 
valid and other modes then the center of mass motions can be excited. The largest oscilla- 
tor strengths are found for the modes closest to the cyclotron frequency. If an external bias 
potential is applied the oscillator strengths of other then the center of mass modes are even 
more enhanced. Crossover phenomena can be observed in this case. For nonequivalent dots, 
even the edge modes have small but non-zero oscillator strengths. 

VIII. ACKNOWLEDGMENTS 

B.P. is an Aspirant and F.M.P. a Research Director of the Flemish Science Foundation 
(FWO-Vlaanderen). This work is partially supported by a NATO-linkage Collaborative 
Research Grant, the FWO-Vlaanderen and the 'Interuniversity Poles of Attraction Program 
- Belgian State, Prime Minister's Office - Federal Office for Scientific, Technical and Cultural 
Affairs'. 



12 



APPENDIX A: SUM RULE FOR THE CLASSICAL OSCILLATOR STRENGTHS 



In this Appendix we present the definition of the classical oscillator strength for a general 
linear mechanical system with weak dissipation and we proof a sum rule for these oscillator 
strengths. We consider a mechanical system interacting with an external force which can be 
described by the following matrix equation 

(j-^-C^X{t) = Af{t). (Al) 

Here symbol X stands for the vector composed of the system variables: coordinates, veloc- 
ities etc. C is the dynamical matrix describing the mechanical system, A is the unit length 
vector (y4^ = 1) characterizing the application of the external force, and the scalar function 
f{t) stands for the external force itself. The power absorption is usually calculated for a 
periodic driving force 

fit) = f{uj) cosM) = Re {/e^-*} . (A2) 
In that case the averaged power absorption is 

PM = -/ dtP{t) = - dtf{t){A^.X{t)) 

= iRe{rH(A^-XM)}, (A3) 

where the coordinate vector Fourier component obeys the following equation 

[iu - C)X{uj) = Af. (A4) 

Note that our set of dynamic equations ( |4.64 |4.6b| ) can be easily reduced to this type of 
equation by just introducing additional variables Y = uX, Z = ujY and constructing the 
resulting vector X = {X, Y, Z} of all those variables. 

Now we use the following dynamical matrix representation 

C = UVU-\ (A5) 

where V is the diagonal matrix composed of the eigenvalues iuji corresponding to the system 
eigenfrequencies. Using that expansion and the substitution u ^ u = uj + irj we obtain the 
solution of equation ([A4|) 



x{uj) = u{io -vy^u-^Af 

= fJ2 U^{^iio - Lu,) - v}-\U-' ■ A), (A6) 



i=l 

and the expression for the power absorption 



1 " 

P{u;) = -fRej:iA^ " ' - vr\uy ■ A) 

"o/^X!? 2' 
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where 



/. = Re{(A^-f/,)(t/ri.A)}, (A8) 

is the oscillator strength. Like in Sec. [V| symbol Ui stands for the i-th vector- column 
of matrix U, and the symbol U^^ is the i-th vector-row of matrix U~^, and the term 
Im |(v4-^ ■ Ui){U~^ ■ ^)} is neglected because in the case of weak dissipation it gives a negli- 
gible eigenfrequency shift. 

The oscillator strengths obey a useful sum rule which follows straightforwardly from 
definition (|A8|), namely, 

n n 
i=l i=l 

= Re(A^ ■ UU-^ ■ A) = (A^ ■A) = l. (A9) 

We obtained a relation which shows that the sum of the oscillator strengths does not depend 
on the properties of the dynamical system. Note we omitted symbol Re in the last equation 
( |X9| ) because the classical equations of type (|A1|) are always formulated with real coefficients. 

The obtained sum rule has a very simple physical meaning. It is just related to the total 
power absorption in the case of the impact external force f{t) = A6{t). Actually, in that 
case calculating the total power absorption 

/oo roo 
dtf{t){A^ ■X{t))= dt5{t){A^ ■ X{t)), (AlO) 
-oo J—oo 

we only need to know the solution in the limit of small time values. That solution can be 
obtained just neglecting the dynamical matrix C in equation (Al) as compared with the 
time derivative. So, the solution reads 

X(t) = A r dt6{t). (All) 

J—oo 

and the total power absorption 

/oo roo rt 1 " 

dtfit)iA^.Xit)) = (A^.A) / dt5it) / dt'5it') = -Y.h. (A12) 
-oo J—oo J—oo Z 

coincides with half of the above defined sum. Thus, the total absorption power does not 
depend on the properties of the dynamical system, and that confirms that the sum of 
oscillator strengths is a universal quantity. 



APPENDIX B: TWO COUPLED OSCILLATORS WITH REPULSIVE 

INTERACTION 

In Sec. VI it was shown that all the spectrum branches show a general behavior, namely 
the frequencies decrease with decreasing interdot distance. Such behavior is caused by the 
repulsive interaction of the electrons belonging to different dots, and can be explained using 
the following simple model. We focus our attention to the center of mass motions of the 
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electrons of the seperate dots. In the case the dots are decoupled (a oo) those motions 
can be treated as two independent oscillators in a magnetic field perpendicular to the plane 
of motion. When the dots are put closer to each other, say, to some distance a ^ Ri, the 
interdot electron interaction can be approximated by 



\2 _i_ r,2 a 2a-^ 



(Bl) 



Here the coordinate corresponds to the lateral center of mass motion of the electrons in 
dot i. The first term in the above expression gives the total energy shift while the second 
one should be inserted into the equations describing the system oscillations. Thus our 
model system consists of two harmonic oscillators with a repulsive harmonic interaction in 
a magnetic field. If also an external clockwise circularly polarized electric field is applied, 
the equations read 



6 6 

mvi + kiTi + -Ti X B r-(ri — = eA^i£^o(cos(a;i)ex — sin(a;i)ey), 

6 

mr2 + /c2ri + -h x B + — (ri - = eA^2-E'o(cos(a;t)ex - sin(a;t)ey), 

C (X 



(B2a) 
(B2b) 



where rii and N2 are the total number of electrons in each dot, and Eq is the amplitude 
of the applied electric field. This linear set of equations can be solved by introducing the 
complex variables Zi = Xi + iyi which leads to the following set of equations 



zi + iUcZi + ulzi — ^q{zi — Z2) 

Z2 + iUJcZ2 + <^2^2 + - Z2) 



eNiEo 
( 

m 

eN2Eo 

( 

m 



-ibjt 



-iwt 



(B3) 



Here Qq — e^/ma^. Now inserting Zi ~ exp(— itut) into the homogenious version of the above 
expressions and mapping the four obtained eigenfrequencies for both polarizations onto the 
positive frequencies we have 



1,2,3,4 



f 00c\ OO1+UJ2 



\\2 J ■ 2 
It is easily seen that the derivative 



— fig ± 



(B4) 



(B5) 



is always negative what confirms that the interdot interaction of electrons causes a negative 
frequency shift, and thus, illustrates the previously obtained negative shift of the spectrum 
branches with decreasing interdot distance. 

The power absorption for this model system can also easily be calculated. It reads 



2m f^ii^- ^i) + T 



(B6) 
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where the oscillator strengths are given by 

It can be checked straightforwardly that the oscillator strengths obey the sum rule Yl'i= 
1. They show the correct qualitative behavior as function of the cyclotron frequency as 
for the coupled dots in Sect. |V1[ 
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FIGURES 



FIG. 1. The layout of the structure. 

FIG. 2. The equilibrium electron density profiles for the system with k = 1 and /ii = /X2 for 
different values of a/Ro- The inset shows the radius of the dots as function of a/Ro. 

FIG. 3. The equilibrium electron density profiles for the system with k, = 2 and fii = H2 for 
the same values of u/Rq as in Fig. ^. The inset in (a) shows the radii of the dots, the inset in (b) 
the number of electrons in dot 1 as function of a/Ro. 

FIG. 4. The frequency spectrum for the system with k = 1. The limiting case a/Ro = oo is 
shown in (a) where the branches are labeled by the corresponding I value. The case a/i?o=0.5 is 
shown in (b). The solid lines correspond to the {a = 1) branches, the dotted lines to the {a = — 1) 
ones. The branches that can be excited in the dipole approximation are indicated by thick curves. 
The corresponding oscillator strengths are shown in (c). 

FIG. 5. The frequency spectrum and the oscillator strength for the system with k = 2 and 
fJ-i = in the limit a/RQ = oo. The notation is the same as in Fig. ^ The correspondence 
between the spectrum branches and the oscillator strengths are indicated by the numbers 1-4. 

FIG. 6. The frequency spectrum and the oscillator strengths for the system with k = 2 and 
^1 = 1^2 foi' the interdot distance a/i?o = 0.5. The numbers give the corresponding frequency 
branches and oscillator strengths. Only 6 oscillator strengths are shown. 

FIG. 7. The frequency spectrum and the oscillator strength for the system with k = 2, 
a/RQ = 0.3 and with an external applied bias potential — fi2 = —0.2. The numbers give 
the corresponding frequency branches and oscillator strengths. Only 7 oscillator strengths are 
shown. 

FIG. 8. The edge modes and the oscillator strengths for the system with n = 2 and /ii = ^2- 
The interdot distance is a/R^ = 0.5. The series ( p.5| ) and ( [4.2| ) are approximated by only the first 
two terms. The numbers give the corresponding frequency branches and oscillator strengths. 
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